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ABSTRACT 

We investigate the reddening and Main Sequence fitted distances to eleven young, 
Galactic open clusters that contain Cepheids. Each cluster contains or is associated 
with at least one Cepheid variable star. Reddening to the clusters is estimated using 
the U-B:B-V colours of the OB stars and the distance modulus to the cluster is 
estimated via B-V:V and V- if .'1^ colour- magnitude diagrams. By main-sequence fitting 
we proceed to calibrate the Cepheid P-L relation and find My=-2.81 xlogP-1.33±0.32 
and Mif =-3.44xlogP-2.20±0.29 and a distance modulus to the LMC of 18.55±0.32 in 
the 1^-band and 18.47±0.29 in the if'-band giving an overall distance modulus to the 
LMC of 18.51±0.3. 

In the case of two important clusters we find that the U-B:B- V diagram in these 
clusters is not well fitted by the standard Main Sequence line. In one case, NGC7790, 
we find that the F stars show a UV excess which if caused by metallicity would imply 
Fe/H~-1.5; this is anomalously low compared to what is expected for young open 
clusters. In a second case, NGC6664, the [/-i?;_S- F diagram shows too red U-B colours 
for the F stars which in this case would imply a higher than solar metallicity. If these 
effects are due to metallicity then it would imply that the Cepheid PL( V) and PL(ir) 
zeropoints depend on metallicity according to gp^jfj ^--^0.66 in the sense that lower 
metallicity Cepheids are intrinsically fainter. Medium- high resolution spectroscopy for 
the main-sequence F stars in these two clusters is needed to determine if metallicity 
really is the ca use or whether some other explanation applies. 



Please see littp://star-www.dur.ac.uk:80^fhoyle/papers.htm] for a version with 



all the figures correctly inserted. 
Key words: Cepheids - open clusters - distance scale - magellanic clouds 



1 INTRODUCTION 

Determining the value of Hubble's constant, Ho, has been a 
challenge to astronomers since the discovery of the univer- 
sal expansion in 1927. It is sometimes argued that we are 
now at the fine tuning stage and many measurements give 
values for Ho which he between the hotly argued values of 
50kms~^Mpc'^ ^ (Sandage) and lOOkms'^Mp c '^ (d e Vau- 
couleurs), e.g. Tanvir, Ferguson and Shanks] (L999) calcu- 
lated Ho=67±7kms^ Mpc~ . However, many of these mea- 
surements are based on secondary indicator methods which 
in turn are dependent on the accuracy of primary indicators 
of distance such as the Cepheid Period-Luminosity (P-L) 
relation. The well-studied LMC P-L relation is usually cal- 
ibrated via the distance modulus to the LMC and the pre- 
viously accepted value was around 18.50. However, this has 



been recently challenged in a paper by Feast & Catchpole 



(1997) who determined the distance modulus to the LMC 
as 18.70±0.1. This small difference in the distance modulus 
causes a 10% decrease in estimates of the Hubble's Con- 
stant. This discrepancy has further motivated us to check 
the Galactic zeropoint of the P-L relation. We do this by 
checking the values of the distance modulus and reddening 
of the 11 Galactic clusters that contain Cepheids via zero 
age main sequence fitting (ZAMS). 

Previous work on measuring the reddening and distance 
to young open clusters which contain Cepheids via ZAMS 
fitting has been done using photoelectric and photographic 
measurements in optical wavebands. It is time consuming to 
observe a large number of stars using photoelectric observa- 
tions as each star has to be observed individually. Photo- 
graphic data can give relatively inaccurate magnitudes and 
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Table 1. Details of the observations. The airmass is the average airmass of the exposure and the exposure time is given in seconds for 
each of the wavebands in column 5. 



colours. However, CCD's now make it possible to observe a 
large number of stars in many different wavebands quickly 
and accurately. Although CCD's h ave alr ea dy been used fo r 
open cluster studie s e.g. Walkei (1985a), Walker (L985b), 
Romeo et al ( 1989| ), these have mainly been carried out in 
BVRI. Recently CCD's with improved 17-band sensitivity 
have become available and [/-band CCD data is included in 
this study. Infra-red imaging detectors are also now available 
and although some of the detectors used here do not cover 
as wide an area as optical CCD's, observing the full extent 
of an open cluster with a mosaic of pointings is a practical 
proposition. 

Until fairly recently, good quality infra-red measure- 
ments of the Cepheids themselves were not available and 
the Cepheid P-L relation has been primarily calibrated in 
the F-band. Laney and Stobie (1993,1994) present infra-red 
along with F-band magnitudes for a large number of South- 
ern Hemisphere Galactic Cepheids. Using data in the litera- 
ture to obtain values for the distance modulus and reddening 
to the clusters they calibrated the Cepheid P-L relation in 
the V and if-band. Any errors in the determination of the 
distance modulus and the reddening in the previous work 
would cause an error in the PL relation as determined by 
Laney and Stobie. 

The layout of this paper is as follows. In section |^ we 
present the observational data and we test the accurac y of 
the photometry and calibration of the data. In section |3.l| 
we describe how the reddenings a nd d istances to the open 
clusters are obtained and in section 3.2 we discuss each clus- 
ter individually. In section ^ we use these values with the 
magnitudes of the Cepheids to calibrate the Cepheid Period- 
Luminosity relation, In section H we discuss the implications 



of the results, particularly for the clusters whose U-B:B-V 
diagrams do not appear to follow the canonical locus. In 
section ^ we draw conclusions. 



2 DATA 

2.1 Observations 

The observations of the Galactic Open clusters were taken 
during five observing runs on the JKT, UKIRT, at CTIO, 
at Calar Alto and on the WHT over a two year period. 
The spread in declination of the clusters and the multi- 
wavelength nature of the study meant that many different 
telescopes were required. 



2.1.1 JKT 

Optical imaging of eight open clusters was obtained dur- 
ing an observing run from the 16/9/1997 to the 22/9/1997. 
The observations were carried out using the 1024 x 1024 Tek- 
tronix CCD with pixel scale of 0.33 arcsec pixel"^. Typical 
seeing was around 1.3". Short exposures of 5s in V, 10s in B 
and 20s in U were observed for calibration purposes but the 
main imaging observations were typically 6xl20s in the V- 
band, 6xl80s in the B-band and 6x300s in the [/-band. Due 
to the Southerly declination of some of the objects, they 
had to be observed at high air mass. However these obser- 
vations were normally used to obtain relative photometry 
and calibration frames were observed at as low an airmass as 
possible or during a later observing run at CTIO. The point- 
ings are given in Table |2[ For calibration purposes, standard 
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Cluster 



Pointing 
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NGC6664 
WZ Sgr 
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NGC6067 
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NGC6823 
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NGC7790 



Star 19 in M adore 
Star 95 in 
Star 5 in 
WZ Sgr 
TW Nor 




Star 136 in Thackeray et al (1962) 
CV Mon 

5" south of TR35 



Star 113 in 



Gluettci 


il99] 


) 


in Aro et all m95n 


1 Romeo et al 


(1989 



Table 2. Approximate pointings for the clusters in the study in 
all wavebands. 



stars from 



Landoll (1992) were used. On the one fully pho- 
tometric night (21/9/97) six Landolt fields were observed at 
regular intervals throughout the night, most of these fields 
containing several standard stars. 



2.1.2 CTIO 

The observations were made using the CTIO 0.9-m dur- 
ing an observing run from the 24/9/98 to 29/9/98. These 
observations were carried out using the 2048x2048 Tek#3 
CCD with pixel scale 0.384 arcsec pixeP^ The Tek#3 CCD 
has low readout noise (4 electrons) and good quantum effi- 
ciency in the [/-band. The average seeing during the obser- 
vations was around 1.4". Approximate pointings are again 
given m Table H Standard E -re gion fiel d s from GrahanJ 
(1982) and standard stars from Landolt (|l992|) were ob- 



served throughout the night. The only photometric night 
was the 28/9/1998 and on this night six E-region standard 
fields and one Landolt standard field containing four stan- 
dard stars were observed throughout the night. Three new 
clusters were observed at CTIO and further observations of 
clusters observed at the JKT were made in cases where the 
clusters had been observed in non-photometric conditions 
only. 



2.1.3 UKIRT 

The infra-red data was mostly observed at UKIRT during 
the four nights 16-19/6/97 using the IRCAM3 near-IR imag- 
ing camera with a 256x256 detector. The pixel scale used was 
0.286 arcsec pixel"^ giving a field of view of 73". To cover a 
sufficient area of each open cluster we therefore had to mo- 
saic images. Generally a mosaic of 9x7 images was observed. 
Each image was overlapped by half in both the x and y di- 
rection so the final image had an approximate area of 6' x 
5'. Observations were generally taken using the ND-STARE 
mode, where the array is reset and read immediately and 
then read again after the exposure which reduces the read- 
out noise to 35e~ . The exposures were 60 x 2 seconds and 
the centre of the mosaic is approximately in the same po- 
sition as the corresponding optical frame. Standards from 
the UKIRT faint standards list were observed throughout 
the nights. Around 10 standards were observed on the three 
photometric nights, some of which were observed early on in 
the night, half way through and at the end of the night. The 



seeing throughout the run was typically O.i 
were also taken in the J and _ff-band. 



2.1.1 C alar Alto 



Observations 



NGC129 lies further north than the declination limit of 
UKIRT so infra red observations were instead taken at Calar 
Alto during another observing run. Observations were done 
using the Rockwell Ikxlk Hawaii detector with pixel scale 
0.396 arcsec pixel"^. This gives a 6. '6 field of view so there 
was no need for the mosaicing technique used at UKIRT. 
The exposure time was 10 x 1.5s, the seeing was b etter than 
1" an d the exposure was centred on star 113 in Arp et al 
( [1959| ). Observations were made in the ifshort-band. UKIRT 



faint standards and standards of Hunt et al (1998) were ob- 
served for calibration. 



2.1.5 WHT 

NGC7790 also lies further north than the UKIRT declina- 
tion limit and so infra red observations were made on the 
4.2m WHT during another observing run. The observations 
were made on the 1/9/1996. The WHIRCAM 256x256 de- 
tector which was situated at the Nasmyth focus and behind 
the MARTINI instrument was used for the observations, 
without MARTINI tip-tih in operation. The WHIRCAM de- 
tector was the IRCAM detector previously used at UKIRT. 
The pixel size was 0.25 arcsec pixel"^ and the field-of-view 
was therefore 64", centred on star E in |R.omeo et al (1989). 
The -Kshort filter was used and UKIRT faint standards were 
observed for calibration. 



All the observations are summarized in Table |l|. The date of 
each observation, the wavebands observed for each cluster 
and the airmass are given in columns 2, 5 and 6 respectively. 
Column 7 gives the exposure time of the frames used for 
imaging. For some of the clusters a calibration frame was 
observed at CTIO and the exposure time of these clusters are 
also given in column 7. Column 8 indicates where a cluster 
was observed on a photometric night and hence where an 
independent zero point was obtained. 

2.2 Data Reduction 

2.2.1 JKT 

Removal of the bias introduced into the data and trimming 
of the frames to remove the overscan region was done on 
all the frames using the IRAF task CCDPROC. At least 
eight 17-band sky flats and six B and F-band sky fiats were 
observed on each of the nights so a separate fiat field was cre- 
ated for every night using a combination of dust and dawn 
sky flats. This was created within FLATCOMBINE using a 
median combining algorithm and a 3cr clipping to remove 
any cosmic rays. The residual gradient in the flat fields is 
around 1%. The task CCDPROC then apphes the fiat fields 
to all the images. The same fiat fields were used in the re- 
duction of the standard star frames. 

Many images of the same cluster were observed. These 
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Cluster 



obs 



^obs 



5„ 



Vohs 



N(stars) Previous Work 



NGC6649 
M25 

NGC6664 
Lynga 6 
NGC6067 
vdBergh 1 
Trumpler 35 
NGC7790 



WZ Sgr 

NGC6823 

NGC129 



0.02±0.025 
-0.01it0.03 
-0.01±0.012 
0.02±0.02 
0.07±0.08 
0.02±0.015 
-0.04±0.02 
-0.009±0.01 



O.OlitO.Ol 
-0.03±0.02 
-0.02±0.01 
0.02±0.02 
0.04±0.03 
0.00±0.02 
-0.01±0.015 
-0.015±0.008 



0.00±0.012 
-0.02±0.02 
-0.01ib0.015 
-0.01±0.015 
-0.06±0.04 
-0.005±0.008 
0.01±0.015 
0.01±0.007 



20 
21 
15 
17 
6* 
24 
15 
22 



O.OOitO.Ol 
O.OOibO.Ol 
0.00±0.03 



O.OOiO.Oll 
0.00±0.005 
-0.005±0.01 



0.00±0.007 
0.00±0.004 
-0.008±0.01 



12 
13 
13 




Table 3. Comparison of the photometry of this study with previous photoelectric data. Note there are 6 stars from Thackeray et al 
( |l962| ) in common with the B and y-band data but only 2 in common with the [/-band data. The last three clusters were observed 
on non-photo metric nights onlv so the previous work was relied upon for calibration. There is a small offset between the zero point of 
NGC129 and [Turner et al (1992) as only the brightest stars were used for calibration. 



were all combined together by aligning the images with lin- 
ear shifts using the task IMSHIFT. Generally these shifts 
were small (a few pixels either way) as the observation were 
done one after each other and in some cases no shifts were 
required. The images were combined using IMCOMBINE 
and were averaged together using a 5(7 clipping. 

2.2.2 CTIO 

The data was obtained at CTIO using the four amplifier 
readout mode. A package called QUADPROC within IRAF 
corrects for the different bias levels in the four quadrants 
of the CCD. The frames were also trimmed usi n g QU AD- 
PROC. In a previous observing run, Croom et a ] ( |1998| ) had 
found that there was a residual gradient of 5% in the dome 
flat fields so sky fiats were used. At least three sky flats were 
observed on each night in the B and l/-band and at least 5 
sky fiats were observed on each night in the t/-band. The re- 
sulting flat flelds were fiat to better then 1%. The equivalent 
version of FLATCOMBINE in the QUAD package was used 
to median combine the flat fields using a 2>a clipping. The 
E-region standards and observations of Landolt standards 
were again reduced in the same manner. 

Multiple images of the same cluster were again com- 
bined using IMCOMBINE with the same settings as for the 
JKT data and where any offset shifts appeared between the 
data frames they were again corrected for using IMSHIFT. 



2.2.3 UKIRT 

The UKIRT data was reduced using a program called 
STRED within the package IRCAMDR. This is a fairly au- 
tomated routine which reads in the data frames, subtracts 
of the dark count and creates a flat field frame by median 
filtering the image frames. Then the program flat flelds the 
dark subtracted object images, corrects for any bad pixels 
and finally creates a mosaic. All the data frames were me- 
dian combined to create a fiat field for each night. There was 
no evidence of a large scale gradient greater than about 1% 
in the fiat fields. 

To create the final image, all the individual frames have 
to be mosaiced together. STRED reads in the offset from 
the data header, however, these offsets were not accurate 



enough. By creating a separate offsets file the mosaicing 
could be done more accurately. To create the offsets file, 
one of the corner frames was fixed and the offset required 
for the neighbour frame were found by eye. This was built 
up over the whole frame, however once one offset had been 
determined, all the other offsets in the x and y direction 
from frame to frame were the same. The offsets for the stan- 
dard stars were more accurate and could be used to create 
the mosaic. 



2.2.i C alar Alto 

Basic IRAF routines such as IMCOMBINE and IMARITH 
were used to reduce the Calar Alto data. A flat field was 
created by median combining all the data frames. A more 
detailed descriptio n of the data reduction can be found in 
McCracken] ( [1999D . We found that the best results were ob- 
tained by first subtracting a sky frame from each image, as 
follows: A sky frame for each individual image was created 
using IMCOMBINE with a 5a clipping to median filter four 
data frames that were local in time to the image frame to 
create a sky frame for each data frame. The sky frame was 
then subtracted off the image before IMCOMBINE was used 
again to combine the data frames, forming one final image 
frame. 



2.2.5 WHT 

The data from the WHT was reduced for us as part of an- 
other project. Dome flats were used to flat fleld the data and 
this was divided into the science frame using IMARITH. Sky 
subtraction was also required. A sky frame was created by 
combining dedicated sky frames observed locally in time to 
the science frame and this was then subtracted from the 
science frame also using IMARITH. 



2.3 Image Alignment 

In order to be able to produce colour-magnitude diagrams, 
the magnitude of each star in all the different wavebands 
is required. To do this, all the frames need to be aligned. 
Aligning the optical frames was easy as there were only lin- 
ear shifts between each waveband. Rather than altering the 
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Figure 1. The airmass coefficient, colour equation and zero point, with the colour equation and airmass correction applied, for the U, 
B and V wavebands for the data from JKT. The open circles in the airmass plot show the standard stars observed at one airmass only. 



data, the alignment was just done by applying small cor- 
rections to the X and y positions of the stars. Aligning the 
optical data with the K-hanA data was more complicated 
though as there were shifts, shears and rotations between 
the frames. We used the IRAF routine GEOMAP to calcu- 
late the best spatial transformation function between any 
two images thus allowing the optical and it-band data to be 
aligned. The mapping was only used to transform coordi- 
nates and we did not perform photometry on the resampled 
images. 



2.4 Photometric Calibration 

2.4.1 JKT 

Observations of the standard stars of Landolt ( 199^ ) were 
taken at regular intervals on each of the six nights at the 
JKT. The CCD frames of the standard stars were reduced in 
the same manner as the data frames with the same flat fields 



etc as discussed in section 2.2.1, The aperture size used to 
measure the magnitude of the standard stars was 15", large 
enough to determine accurately the total magnitude of the 
star but not so large that sky subtraction errors dominate 
the magnitude measurement. 

Only the sixth night (21/9/97) was fully photometric 
and the zero point, airmass coefficient and colour equation 
are shown in Figure ^ and given below with the rms scatter. 

f/jkt = (/idt + 3.76 + Q.4Qsec{z) - 0.063((/idt - -Bidt) ± 0.038 

Bjkt = Bidt + 2.00 + 0.27sec(z) - 0.013(Bidt - Mdt) ± 0.033 

Vjkt = '^dt + 1.99 + 0.19sec(2) ± 0.026 

where the subscript Idt stands for the Landolt standard star 



magnit ide and the subscript jkt stands for the instrumen- 



tal magnitude, z represents the zenith distance. The errors 
on the airmass are ±0.0025, 0.0012 and 0.0008 in the U, B 
and V and the errors on the colour equation are ±0.0007 
and 0.0005 in the U and B-bands respectively. There are 
no Landolt magnitudes available for the data frames so the 
colour term has to be translated in instrumental magnitudes. 
The colour term is negligible in the F-band calibration so 
the instrumental V-band magnitudes come directly from the 



above equations. The 5idt- Vidt and (/idt-Bidt colours are 
given by 

(Bidt - Vidt) = 1.013((Bjkt - Vjkt) - 0.01 - 0.08sec(2)) 

(f/idt - Bidt) = 1.067((f/jkt - Bjkt) - 1.76 - 0.22sec(2)) 

we assume that the contribution from the colour term in 
the calibration of the B-band is negligible when determining 
the I7idt-Sidt colour. These exact colour terms are used to 
correct the instrumental magnitudes in order to make the 
colour-colour and colour-magnitude diagrams in Figures |ll], 
|l| and |l|. 

The clusters NGC7790, NGC6664 and Trumpler 35 
were observed on the one photometric night. Short exposure 
observations of NGC6649 and M25 were made at CTIO in 
order to obtain an independent zero point for these frames. 
Around 20 bright (brighter than V~15), fairly uncrowded, 
unsaturated stars were taken as standard stars to identify 
the relation between the zero point from the JKT data and 
from CTIO data to an accuracy of 0.01 mags. The agree- 
ment between the zero point found via this method and 
previous, photoelectric calibrations is good (see Table ^ 
with only small offsets in each case. The remaining three 
clusters were observed in non-photometric conditions only 
so previous work had to be relied upon for the calibration. 
Fo r NGC6 8 23 th e photoelectric observations from Table 1 
of Guettei (1991) were used. Some of these stars were sat- 
urated on the CCD frame and the area of overlap between 
the two images was not identica l but th i rteen stars were 
suitable for calibration purposes. Guetter (1991) compares 



his photoelectric data with that of previous work and finds 
good agreement. Two sources of photoelectri c data are avail- 
a ble for the cluster NGC129, ( |Arp et ajp^Hg] ) and ( [Turner et 



a] 1992 ). There are 9 stars in common with the Arp photom- 



etry and 13 stars in common with the Turner photometry. 
For these samples of stars, we find that the t^-band zero 
point obtained from Arp is 0.04±0.03 mags brighter than 
that of Turner. In the B-band the difference is 0.02±0.01 
mags in the sense that Arp is brighter than Turner and 
there is 0.01±0.01 mag difference in the same sense in the 
y-band. These differences are mainly caused by the stars in 
the sample with F fainter than 14 mag. We therefore use the 
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Figure 2. The airmass coefficient, colour equation and zero point, with the colour equation and airmass correction applied, for the U, 
B and 1/ wavebands for the data from CTIO. The open circles in the airmass plot show the standard stars observed at one airmass only. 



average of the brightest two stars from Arp and Turner to 
calibrate NGC129. The zero point from this method agrees 
very weU with the zero point obtained using Turner's pho- 
tometry, which is shown in Figure |l6[ FinaUy the photoelec- 
tric work of Turner (1984) was used to calibrate the cluster 
containing the Cepheid WZ Sgr. 

The magnitudes of the stars on the data frames on all 
nights were measured using a 5" aperture. The standard 
star magnitudes were measured using a 15" aperture so for 
the photometric night (21/9/97) an aperture correction had 
to be applied to the data. The correction in the U, B and 
y-bands are -0.165±0.02, -0.11±0.025 and -0.11±0.024 mag- 
nitudes respectively as determined by comparing the mag- 
nitudes of the standard stars with a 5" and 15" aperture. 
No aperture correction was required for those cases where 
the calibration was tied to a photometric sequence in the 
cluster itself. For the clusters NGC6649 and M25 the CTIO 
aperture correction was required. 



24.2 CTIO 

E-region standards from Glraham ( 1982| ) and Landolt stan- 
dards were observe d for the phot o metr ic calibration of the 
optical CTIO data. VIenzies et a\ ( 1991 ) compared the zero 
points and colour differences found from using the two dif- 
ferent standard star studies and found the offsets between 
the two to be small, 0.004±0.0095 offset in the sense E re- 
gions - Landolt. and similar sized offsets in the U-B and B- V 
colours. Any offsets are within the quoted error. As before, 
these standard frames were reduced in the same manner as 
the data frames and a 15" aperture was used to determine 
the magnitude. Again, only one night was photometric and 
this was the last night (28/9/98). The zero points, airmass 
coefficients and colour equations for each waveband, [/, B 
and V, are shown in Figure ^ and given below again with 
the rms scatter. 



where the subscript std stands for the E-field or Landolt 
standard magnitude and ctio stands for the instrumental 
magnitude. The errors on the airmass are ±0.003, 0.0011 and 
0.0009 in the U, B and Fand the errors on the colour equa- 
tion are ±0.0006, 0.0007 and 0.001 in the U, B and F-bands 
respectively. There is generally good agreement between the 
values for the airm ass coefficients and colour equations found 
in this work and in Croom et a] ( 1995 ) . Again, a 5" aperture 



f/ctio = U,td ± 4.61 + 0.47sec(^) - 0.036([/std - Bstd) ± 0.035 
Bctio = Bstd ± 3.15 ± Q.2Uec{z) ± 0.099(B,td - K^td) ± 0.015 
Ktio = Vstd + 2.93 ± Q.l2sec{z) - 0.018(Bstd - Ktd) ± 0.007 



was used for the data frames so an aperture correction for 
the photometric night was required. The aperture correc- 
tions, in the U, B and F-bands, are -0.17±0.03, -0.17±0.02 
and -0.13±0.02 magnitudes respectively. Again, for the data 
frames the colour terms have to be found in terms of CCD 
magnitudes rather than standard magnitudes. The colour 
term in the B-band is in this case non-negligible so F-band 
magnitudes have to be used in the [/-band calibration. The 
standard colours are given by 

(Bstd - Ktd) = 1.088((Bctio - Ktio) - 0.22 - 0.09sec(2)) 

{Um - Bm) = 1.037(((7ctio - Bctio) - 1.46 - 0.26sec(z) 
+0.099(1.088((Bctio - Ktio) - 0.22 - 0.09sec(2)))) 

again these colours are used in all to calculate the instru- 
mental magnitudes which are used in the colour-colour and 
colour-magnitude diagrams in Figures ^ and ^| 

This night provided independent zero points for the 
clusters NGC6067, Lynga 6 and vdBerghl, which were 
not observed at JKT, and also provided a zero point for 
NGC6649 and M25, which were only observed in non- 
photometric conditions at the JKT. 



24.3 UKIRT 

The standards stars observed at UKIRT were taken from 
the faint standards list available from the UKIRT Web page. 
The standard stars were observed as a mosaic of five frames, 
a central frame with an overlapping frame in each direction. 
A 15" aperture was used to determine the magnitude of the 
standard. Three of the nights were photometric, nightl,2 
and 4 (16,17,19/6/97) so all the clusters were observed in 
photometric conditions. 

The zero points for each of the nights are shown in 
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Figure 3. The i<'-band zero points for tlie three photometric 
nights at UKIRT. The top panel is for night 1, (16/8/97, fiUed 
circles), centre for night 2 (17/8/97, open circles) and the lower 
panel for night 4 (19/8/97, filled triangles). 



Figure 4. The airmass correction and colour equation for the 
if-band data taken at UKIRT. The symbols are the same as in 
Figure H. 



Figure q and the airmass coefficient and colour equation are 
shown in Figure ^ and given below for each night. 
16/06/1997; 

Kv = ifstd + 6.94 - 0.082sec(z) + 0.005(J,td - K,td) ± 0.024 
17/06/1997; 

Kv = ifstd + 6.95 - 0.082sec(2) + 0.005(Jstd - K,td) ± 0.033 
18/06/1997; 

Ku = ifstd + 6.96 - 0.082sec(z) + 0.005(J,td - K,td) ± 0.037 

where again the subscript std refers to the standard stars. 
The difference between each of the cahbrations is just a smaU 
change in the zero point. The error on the airmass coeffi- 
cient is ±0.002 and the error on the colour term is ±0.0003. 
The airmass coe f ficien t agrees weU with the values given in 
Krisciunas et al (1987) and those found on the UKIRT web 
page. An aperture correction of -0.19±0.018 magnitudes is 
required for the magnitudes of the data frames measured us- 
ing a 5" aperture on each night. As the colour term is very 
small, it was assumed negligible in the V-K:V CMD's. 
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2.4.4 Calar Alto and WHT 



UKIRT faint standards and Hunt et a] ( 



199S ) standards were 



observed in order to calibrate the Calar Alto and WHT data. 
The calibration was provided for us by Nigel Metcalfe as part 
of another project (see McCracken, Metcalfe and Shanks, in 
preparation, for more details). 



Figure 5. An indication of the accuracy of the f7-band photome- 
try based on results for simulated stars (see text for details). The 
large errorbar indicates the total error (Poisson errors, crowding 
errors and also read noise), the smaller, wider erro rbar indicates 
the Poisson error found as described in section |2.5[ paragraph 2. 
The total error at 17th magnitude is only ±0.036 magnitudes and 
down to 18th magnitude the error is less than ±0.1 magnitudes. 



2.5 Photometry 

Automated aperture photometry was done using PHOT 
within IRAF's DAOPHOT package. A smaU, 5", aperture 
was used to minimise any crowding problems. PHOT was 
also used to obtain the magnitudes of the standard stars, 
the magnitude of each standard star was measured individ- 
ually using a 15" aperture. 

First we establish the depth of the photometry. We 
define the limiting depth of the observations to be where 



the Poisson error in the electron counts is 5%, i.e. when 
y/Nobj ± N sky /Nob] = 0.05. The depth of the JKT data is 
19.9 mags for a 1200s (/-band exposure, 19.7 mags for a 360s 
B-band exposure and 20.4 mags for a 240s F-band exposure. 
The depths of the CTIO data are similar. The depth of a 
typical i^'-band exposure of 120s observed was around 19.2 
magnitudes. 

The accuracy of the zero points were tested by com- 
paring our photometry to photoelectric observations in the 
literature. Table shows the residual when our magnitudes 
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for typ icall y 10 stars from the calibrated data frames (see 
section 2.4) are compared to the magnitudes found from pre- 



vious photoelectric studies. In all cases except for NGC6067, 
the resi dual is small (less than 0.03 mags). The errors in the 



comparison of our magnitudes to other photoelectric ob- 



servations are slightly smaller than the errors in the zero 
point found from the standard stars in this work. This is 
probably because the stars considered when comparing the 
previous work were slightly brighter. Unfortunately for the 
cluster NGC6067 there is little photoelectric data available 
and only 6 stars can be compared in the B and F-bands 
due to saturation and different parts of the cluster being 
observed. In the [/-band there are only 2 stars in common. 
Given that the offsets between the work here and previous 
work are small for all the other clusters, we assume that the 
zero point obtained for NGC6067 is accurate. 

Finally, we test how the crowding of the field effects 
the photometry. Twenty "simulated stars" of each magni- 
tude shown in Figure ^ were placed in a 300s 17-band image 
of the cluster M25 then PHOT was used to determine how 
well the magnitudes could be recovered. Shown in Figure 
^ are the results. The x axis shows the true magnitude of 
the stars and the y axis shows the deviation of the mean 
magnitude from the true magnitude. This plot shows the 
total error (Poisson errors, read noise and also errors due to 
crowding) at each magnitude (full errorbar) and the Poisson 
error, found as described above, as the smaller, wider, error- 
bar. The total error for this exposure of 300s is less than 
0.03 magnitudes down to U=17 mags. For a [/-band expo- 
sure of 1800s, assuming that the crowding errors remain the 
same and then correcting the error shown in Figure ^ for 
the reduced Poisson error, the total error is estimated to be 
±0.17 at (7=20, reducing to ±0.08 at 18th mag and ±0.03 
at 17th mag. Figure ^ also shows that although the errors 
increase for fainter magnitudes there is no systematic trend. 



3 REDDENING AND DISTANCE 
3.1 Method 

In order to work out the Cepheid Period-Luminosity (P-L) 
relation, the distance to the cluster needs to be known. How- 
ever, there is generally significant dust absorption along the 
line of sight to the cluster which must be corrected for. The 
method for determining both of these parameters is done 
via ZAMS fitting to colour-magnitude and colour-colour di- 
agrams. Th e ZAMS used in this study is a combination of 
AUerj ( 197^ ) for the optical data, the intrinsic colours of near 



IR-band stars from the UKIRT Web page (Tokunaga 



199J 



for the K-hand data and nearby local stars taken from the 
Strasbourg Catalogue to give an indicatio n of the accept - 
able spread in t he Z AMS. The ZAMS of [Turnei] ( |l979b| ) 
and Mermillio^ (1981) have also been tested and would give 
the same results as the ZAMS of lAUeiJ (|1973|). 



The reddening is obtained using the I7-5/5- V diagram 
which is independent of the distance to the cluster. The U- 
B:B-V diagram has always been the favoured method for 
estimating the reddening; however, previously, the accuracy 
of the reddening determination was limited by the depth 
of the [/-band photoelectric photometry. The improved U 
sensitivity of the current generation of CCD's should allow 



a potential improvement in the accuracy of the reddening 
estimated from U-B:B- V diagrams and this is the route we 
have adopted here. 



The r eddening law assumed in this work is from sharp- 

les3 (|1963|). 



EjU-B) 
E(B - V) 



= 0.72±0.05E(B- V) 



(1) 



We choose to fit the ZAMS to the ridge-line of the O 
and B stars rather than the least reddened envelope, be- 
cause it helps take account of differential reddening in some 
of the clusters. There is evidence for differential reddening 
in NGC6823 and TR35 (see Figure ^l]) as the main sequence 
in the colour-colour diagram is substantially broadened. By 
fitting to the centre of the data we measure the average red- 
dening for the cluster which we can then apply to colour- 
magnitude diagrams which are uncorrected for differential 
absorption to obtain distances. To determine the error we 
measure the standard deviation of the O and B type stars 
from the ZAMS via least squares fitting. The errors quoted 
on the values for the reddening (see Table ^) are typically 
0.1 mags and include any error in the calibration. In quite 
a few cases, the ZAMS does not fit the colour-colour data 
well over the whole range of B-V colours. This is particularly 
problematic in some cases and these cases are discussed be- 
low. There is also the problem that the Cepheid could have a 
different reddening to the cluster, caused by differential red- 
dening across the cluster or by the location of the Cepheid 
away from the cluster. This is discussed further in section ^ 

We use both the B- V: V and V-K: V colour- magnitude 
diagrams to determine the cluster distance. F-ii"; V diagrams 
have the advantage that the slope of the ZAMS is flatter 
than at B- V: V, possibly allowing more accurate distance es- 
timates but V-K: V diagrams are available for only 7 of the 
clusters and the scatter in these diagrams is greater than in 
the B-V:V CMD's, particularly in the case of NGC6664 as 
there were difficulties aligning the V and K frames due to a 
shift in the telescope position half way through creating the 
if-band mosaic for this cluster. There are also few points on 
the V-K:V diagram for NGC7790 due to the small size of 
the IRCAM detector. This diagram was also made for us be- 
fore the F-band observations for this cluster were available 
as a test for the feasib ility of this p r oject . The IZ-band data 
therefore comes from [Romeo et al (1989), however there is 
a good match between our l/-band data and that of Romeo 
et al. (see Figure The distance modulus which best fits 
the B-V:V CMD around the position of the AOV stars, us- 
ing the method of least squares, is taken to be the distance 
modulus, /io, of the cluster. The errors on the distance mod- 
ulus were found by measuring the standard deviation away 
from the ZAMS of AOV type stars over the range -Q.l<^B-V 
<0.1 in the dereddened ZAMS which covers a range of ap- 
proximately 4 mags in V. The distance modulus found from 
the B- V: V CMD is then checked against the V-K: V CMD 
for the 7 clusters with such a CMD for consistency. In all 
cases the distance modulus found from the B- V: V diagram 
was consistent with the V-K: V CMD within the errors. 

No attempt has been made to remove foreground and 
background stars. Only stars which lie clearly off the main 
sequence (off in B- V by more than 1 mag for example) were 
removed. There is no clear recipe for how to remove the 
contaminating stars from the colour-magnitude and colour- 
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Cluster 



E(i3-V0ciu 



Mo 



NGC6649 
M25 

NGC6664 
WZ Sgr 
Lynga 6 
NGC6067 
vdBerghl 
Trumpler 35 
NGC6823 
NGC129 
NGC7790 



1.37±0.07 
0.49±0.08 
0.66±0.08 
0.56±0.20 
1.36±0.17 
0.42±0.02 
0.90±0.18 
1.19±0.10 
0.85±0.09 
0.57±0.06 
0.59±0.05 



1.35 
0.48 
0.64 
0.57 
1.34 
0.35 
0.77 
0.92 
0.75 
0.53* 
0.64* 



11.22±0.32 
9.05±0.43 
11.01±0.37 
11.15±0.49 
11.1±0.45 
11.17±0.35 
11.4±0.65 
11.3±0.53 
11.20±0.55 
10.90±0.37 
12.70±0.43 



11.28 
9.03 
10.40 
11.22 
11.43 
11.13 
11.36 
11.56 
11.79 
11.24* 
12.39* 



Table 4. The values in the abov e Table come from the work here 
and from Lancv fc Sto bic (1993). The values marked * come from 



Feast & Walker 



(1987) as the clusters NGC129 and NGC7790 are 



not included in the studies of Laney and Stobie. 

colour diagrams so the data shown in Figures ^ and 
[l^ contain non-cluster members. The observations in this 
study only cover a field-of-view of ~ 6' and are pointed at 
the cluster centre so contamination may not be as much of 
an issue as for the wider field photographic plates. As an 
example, we discuss this further for the cluster NGC7790 in 
section p. 4 

All of the U-B:B- V diagramMhe B- V: V and the V-K: V 
diagrams are given in Figures |l^, |l2| and |l3| respectively 
and our results for the reddening and distances obtained 
are given in Table ^ togethe r with previous results as sum- 
marised by Laney & Stobie (1994). All of the clusters are 



individually discussed below. 



3.2 Discussion of Individual Clusters 



NGC6649 has been studied prev iously by. for example . 
Madore fc van den Bergh| (|l975| ) and [Walker fc Laneyj ( ^986| ) . 



The agreement between the photometry of this study a nd 
the photoelectric data of Madore & van den Bergh (1975) is 
good (see Table H) , with only small off sets in the U and B- 



band. [Madore fc van den Bergh| (|1975|) find E(_B-V) = 1.37 
(no quoted error) for the r eddening towards the clust e r. The 
distance modulus /!□ from Madore & van den Bergh (1971:) 
is 11.15±0.7. 



Tu rner ( 1981 ) uses th e photometry of Madore & van den 



Bergh) ( |l975D and that of [Talbert| ( |l975| ) to study NGC6649 
and find the cluster suffers from differential reddening. How- 
ever for stars close to the cluster centre a value of E(_B- 
V)=1.38 is appropriate. The distance modulus is found to 
be 11.06=fc0.03 when i n divid ual stars are dereddened. 

Walker & Laney ([l986[) used U, B and F-band CCD 



data to stu dy NGC6649. Agreement b e tween the photoelec- 
tric data of [Madore & van den Bergh (1975) and Walker & 



Laney ( 198( ) was found to be b etter than 0.03 mags in the 
V-band. Walker & Laney (1986) do not measure the redden- 



ing of the c luster due to the claims of differential reddening 
by Turner. Walker & Laney (1986) deredden each star indi- 



vidually to find a distance modulus of 11.00±0.15. 

Rather than correct for differential reddening, we fit 
the ZAMS hue to the centre of the U-B:B- V and B- V: V di- 
agrams to try and measure the average values of the redden- 
ing. We find E(5-T/)=1.37±0.07 and /Xo=11.22±0.32, both 
values are consistent with previous work. The distance mod- 



ulus used in [Laney fc Stobie| ( |l994| ) is slightly higher, 11.278 
and the reddening, E(_B- V), slightly lower, 1.35 but again 
these values are well within the errors. 



Barrell (1980) found the radial velocity of the Cepheid 
V367 Set to be -20±6 kms~^ and the radial velocity of the 
cluster NGC6649 to be -14±5 kms"^ This is taken to be 
evidence for the cluster membership of the Cepheid. 

M2 5 The photometry used in this study and that of 
Sandage ( |L960| ) is compared in Figure [14 and ha. The agree- 
ment in all the wavebands is good, less that 0.03 mags dif- 
ferent from that of Sandage. 

M25 has a U-B:B- V diagram where the ZAMS fit is 
good over a wide range of B-V colours. The value for E(_B- 
V ) of 0.49=fc0.08 is in excellent agreement with the work 
of 



Sandage (1960) who obtained 0.49±0.05. Johnson 



stated that E(_B-V) lay in the range 0.4 to 0.56 and [van de: 



Bergh] (|l97^ obtained E(B- VO=0.51±0.01. 

The distance modulus obtained for this cluster is 
9.05±0.43 which is slightly higher than that of Sandage, who 
found /io=8.78±0.15. The difference in the distance modu- 
lus is probably due to where to fit was made, we fit to the 
A type stars to obtain 9.05±0.43, Sandage's work contains 
brighter stars which may lie slightly off the main sequence. 



Wamp ler et al (1960) obtained 9.08±0.2 and van den Bergh 
( 1978| ) found ^0= 9.0±0.3 which agree well. There is also 



good agreement between our values and those used by Laney 
fc Stobiej ( [l994|) , as given in Table ^. 



"FSast ( 1957 ) studied M25 by measuring radial velocities 
and spectra for stars within M25. He obtained radial veloc- 
ities of around 4 kms~^ for the Cepheid U Sgr and for 35 
stars in the cluster M25 indicating that U Sgr is a member 
of M25. 

NGC6664 has been studied previously by 



(1958) 



The agreement between his photometry and ours is good, 
with only small offsets between the two data sets as given 
in Table |[ 

Figure ^ shows that the cluster NGC6664 has an 
anomalous I7--B;-B- F diagram. This is a clear case where we 
fit the O and B type stars rather than trying to fit the F and 
G type stars as the F and G type stars maybe more affected 
by metaUicity. The value for E{B- V) is then 0.66± 0.08 which 
is slightly larger than 0.6 obtained by Arp (1958). No error 
is quoted by Arp. In the previous work by Arp, the obser- 
vations were not deep enough to see if the anomalous shape 
of the (7-_B;i?- 1/ diagram would have been detected or not. 
Unfortunately the only previous source of photometry for 
more than a handful of stars is that of Arp so no other zero 
point comparisons can be made. 

The distance modulus obtained is 11.01±0.37 which is 



larger than 10.8 from Arp (1958), mostly due to the in- 
creased value for the reddening estimated here. 

There is an larger dis crepancy betwee r i the value of the 
distance modulus used by Laney fc Stobic ( 1994 ) who quote 
a reddening of 0.64 but a distance modulus of 10. 405. It is 
not clear why their distance modulus is so low as Feast fc 
Walkeil (|1987) use 10.88. 



The radial velocity work of Kraft ( 1958 ) shows EV Set 
is a member of NGC6664. 

WZ Sgr The first problem with WZ Sgr is t hat its 
memb ership of an open cluster is questionable. Turner 



(1984) discusses the membership of WZ Sgr to the cluster 
C1814-190 in some detail and concludes that the strongest 
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evidence for membership of WZ Sgr to an open cluster comes 
from the fact that often the Cepheid in a cluster is around 
4 magnitudes more luminous than the B-type stars on the 
main sequence. This essentially means that the a ge of the 
Cepheid is co nsistent with the age of the cluster. Feast & 
Walkei{|(|l987|) assume WZ Sgr to be a cluster member. 

The cluster C1814-190 is also only comparatively 
sparsely po pulated with o nly around 35 members brighter 
than B-- 



-16 



Turner (1984). There is also patchiness in the 



dust obscura tion which would cause differential reddening 
(Turner 1984). These two factors could possibly explain the 



odd shape of the [/--B/B- F colour-colour diagram particu- 
larly in the range 0.5 < B-V < 1. Contamination from fore- 
ground and background stars could also be the source of 
the unusual [/-_B;_B- 1/ diagram. By measuring the reddening 
of the O and B type stars though we get a value fo r E(_B- 
'|/)=0.56±0.20 which is in agreement wit h iTurneil (|l984|) 
and c onsistent with E(_B-V)=0.57 used by Laney & Stobie 

(Fggi). 



The B- F; 1/ diagram is surprisingly tight, giving a dis- 
tance niodulusof 11. 15±0. 49 which is again in agreement 
with [Turncil ( |l984|) who obtained 1L16±0.1. The distance 
modulus used by Laney fc Stobie (1994) is slightly larger, 



11.219 but again well within the errors of our result. 

Lynga 6 The photometry o f Lyn ga 6 has been checked 
Vogt ( 1975) (shown in Figure [l4 



Moffat 



against that of 

a nd |l5|) and also against that of van den Bergh & Harri: 
971 



Offsets between the different sets of photometry are 



less than 0.03 mags in each waveband. 

Like NGC6649, Lynga 6 is also very heavily reddened 
which makes it very difficult to observe stars over a large 
range of B-V and U-B colours. E(_B-V) is estimated to be 
1.3 6±0.17 which is consistent wit h previous measurements 
by v&n den Bergh fc Harris (1976) who obtained 1.34±0.01, 



by Madore ( 1975) who obtained 1.3 7±0.03 and the value of 

I. 34 used by |Laiiey fc Stobitj ( [l994[ ). 

The distance modulus obtained here is around 

II. 10±0.45 which is consistent with 11.15±0.3 obtained 



by Walkei ( 1985a ) but is discrepant with the value used 
by Laney & Stobie (1994). They have a distance modulus 



of 11.429 whi ch agrees reasonably well with the value of 
Madord (|l975|) who obtained (V'-Mv)=16.2±0.5 for the ap- 



parent distance modulus which, if a value of 3.2 is assumed 
for the extinction coefficient, roughly implies 11.8±0.5 for 
the distance modulus (see section 4| for more details of the 
extinction coefficient used here). The difference between the 
value of Madore and the value for the distance found here 
is that Madore tended to fit the edge of the ZAMS. 

The Cepheid TW Nor lies close to the centre of the 
cluster Lynga 6 and has a very similar value of the reddening. 
This is taken a s evidence of t he membership of the Cepheid 
to the cluster ( Walker 1985a). 

NGC6067 There is very little photoelectric data for 
this clu ster. The B and F - band s have been compared to the 
data of Thackeray et a] (1962) but there are only 7 stars 



in common and in the (7-band there are only two stars in 
common. The comparison shows that the zero point used 
here is at least consistent with previous the previous work. 
As the photometry obtained at CTIO for other clusters such 
as Lynga 6 and vdBerghl agrees well with previous results, 
we have to assume that t he photometry for NGC6067 is 
also good. Walker (1985b) finds good agreement between 



th eir B and F-band CC D data and the photoelectric data 
of [Thackeray eTai] (|l962|). 



NGC6067 has the lowest value for the reddening, of the 
clusters in this study, with E(5-l/) estimated as 0.42±0.02. 
The (7--B;i?- 1/ diagram presented here has a main sequence 
which agrees fairly well with the ZAMS, although there is a 
spread around B- V= 0.8. The value of 0.42 i s sligh tly higher 
than previous values, Goulson and Caldwcl^ (1984) obtained 



Thackeray et a| (1962) obtained 0.33 with no 



0.35±0.1 and 
quoted error. 

As the reddening has been estimated to be slightly 
higher than previously, the distance modulus is also in- 
creased t o 11.17±0.35 as o ppose d to 11.05±0.1 from Walker. 
Ifowever Thackeray et a^ (1962) found the distance modulus 
of the cluster to be around 11.3. This estimate is higher than 
our estimate, despite a smaller measured reddening, as the 
ZAMS fit was made t o the edge of the CMD. 

Laney & Stobie (1994) use 0.35 for the reddening and 
11.13 for the distance modulus to the cluster. 

The membershi p of the Cephe ids to this cluster is dis- 
cussed in detail by Eggen (1983). Eggen noted that the 



Cepheid V340 Nor is centrally located in the cluster and 
has the same reddening as the cluster so is assumed to be 
a member. The Cepheid QZ Nor lies a t a distance of two 
cluster radii out from the cluster centre ( Walker 1985b) but 
i s still assumed to be a cluster member by Feast & Walker 
( [1987D 

vdBerghl Ou r pho t omet ry is tested against the pho- 
toelectric data of of Arp (196C) in Table ^. There is no sig- 
nificant offset between the t wo data sets a nd n o evidence 
of any scale dependent error. Turner et a yi998D has com- 
pared their photometry to that of Arp (1960) and find good 
agreement. 

The cluster vdBerghl was only given a short exposure 
of 300s in U and 90s in B and V. The U-B:B- V diagram 
is therefore not very well populated at faint U magnitudes. 
There is a fair amount of scatter in the [/--B/B- V diagram 
for this cluster, the average value is E(_B- 'l/)=0.9 ±0.18. This 
is l arger than the value of previous estimates 
( [1998D obtained a minimum value of E( i ?- V) 



Turner et 



=0.66, this 



clearly fits the the edge of the B-type stars. Arp (1960) ob- 



tained 0.76 (no error) but the spread in the U-B:B-V dia- 
gram is such that the larger value would also have b een ac- 



ceptable. 0.77 w as used for the cluster reddening by Laney 
k Stobie! (|l994| ) 

The distance modulus, /io obtained here is ll.4zfc 0.65. 



Turner et a] 



This is 

mostly due to the measuremen t of in creased reddening. The 



3 ( |1998| ) 



This i s larger than previous re sults, 10.94 o b taine d by Arp 
( [1960| ) and 11.08 obtained by 



value used by Laney fc Stobie (1994) was 11.356 



The me mbership of CV Mon to the cluster has been de- 
termined by Turner et al (1998) using radial velocity mea- 



surements, evolutionary arguments and by it's location in 
the cluster. 

Trumpler 35 There is very little photoelectric data 
available for this cluster. We compare our phot o metry 
against the photoelectric observations of Hoag et al (1961) 



and find reasonable agreement in the B and F-bands. The 
agreement in the ?7-band is less good but as there is only 
one source of comparison and the JKT photometry appears 
to agree well for other clusters we suggest that the [/-band 
data photometry is accurate. 
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Th e data in the t/-_B.-j?- V diagram for the cluster Trum- for the red dening and the distance modulus use d by [Feast 



pier 35 (TR35) shows quite a large spread. This is probably fc Walker (1987) agree better with the values of 



caused oy differential reddening. Rather than trying to cor- 
rect for this, a mean value of E(_B- V)=1.19±0.10 is taken for 
the reddening. Shown on the TR35 panel in Figure |l^ are 
two l ines. Th e dashed shows the value of E(_B- 1^=1.03 taken 
from Turner (1980) and the solid line shows the value of 1.19 



Turner et 



adopted here. Using this, the distance modulus to the clus- 
ter is then estimated as 11.3±0.53 by fitt ing up 1:he cer itre of 
the data. This is different to the value of Turner (198C) who 
obtained 11.6±0.16 as Turner fitted the edge of the B-V:V 
diagram rather than the c e ntre. The values for the reddening 
used by [Laney fc Stobi^ ( |l994| ) (E(B -V)=0.92, /Xo = 11.56) 



are in better agreement with those of Turner ( 1980 ) rather 
than those obtained here. 

The membership of RU Set to Trumpler 35 using the 
reddening and evolutiona ry status of t he Cepheid is dis- 
cussed and supported by Turner] (1980|). H owever RU Set 
does he 15' away from TR35 ( |Tumer| [l980|) . 

NGC6823 was only observed in non-p hotometr i c con - 
ditions. We use the pho toelectric data of Suetter (1991) 



for calibration purposes. Guetter (1991) has compared his 



photometry with that of Hiltner (1956) and Hoag (1961) 
and finds that there are only small offsets of les s than 
0.04ma gs in U-B between the different data sets. Turnci 
( 1979a ) finds good agreement with the photometry of Hilt- 
i ier (1 956) which agrees well with the photometry of Guetter 
(1991) used for calibration here. 



NGC6823 suffers from differential reddening (Turner 

|l979a| ), perhaps to an even greater extent than Trumpler 35. 
Again shown in the panel for the cluster NGC6823 i n Figure 
|ll| are t wo lin es. One is for F,{B-V)—0.53 taken from Feast & 
Walkei[l( |l987| ) and the other is E(5-l/)=0.85 which fits the 
centre of the B-type stars. The value of E(_B- V)=0.85±0.09 
is assumed here. 

The distance modulus with this reddening i s then 
1L20±0.55, lower than the value of 11.81 found by Turner 



(1979a) who again fi tted the edge rat h er th an the centre of 



the B- y.- y diagram. Laney fc Stobie (1994) assume a red- 



dening of E(_B-l/)=0.44 and = 11/787. These values again 
agree fairly well with those of Turner ( 1979a ) rather than 
the values found here. 

The Cepheid SV Vul is only thought to be associated 
with the cluster NGC6823 ( [Feast fc Walker) |l987i) as it lies 
a few arcmins away from the cluster centre. 

NGC129 was only observed in non-photometric condi- 
tions so previous work had to be relied upon for the calibra- 
tion, as discussed in detail in section 2.4.1. The reddening 



for this cluster is E(g- V)=0. 57±0.06. This value is slightly 
larger than previous values. Turner et al (1992) obtained 



0.47 but fits to the least reddened edge of the O and B star s 
rather than to the centre of these stars. Arp et a 

] ( |l95£| ) 

found 0.53 (with no error) for the reddening which is in 
agreement with the value found here. 

The distance modulus obtained assuming a value of 0.57 
for the reddening is /^o = 10.90±0.37. This value is lower than 



the value of 11.11 obtained by Turner et a (1992) who again 



fits the e d ge of the main sequence rather than the centre. 
Arp et a] ( |l95S| ) find 11.0±0.15 for the distance modulus of 
the cluster, in agreement with t he value found he r e. PL Cas 
is not included in the study by Laney & Stobie (1994) due 



alj ( 199S )_rathcr than those obtained here. 

Kraft (1958) found a measurement of -14±3 kms^^ for 
the radial velocity of the cluster NGC129. He found that 
the Cepheid itself had a radial velocity of -llkms"^. Given 
that the error on any individual measurement is estimated 
to be around 1.5kms~^, DL Cas is assumed to be a member 
of NGC129. 

NGC7790 The 1/ diagram for this cluster ap- 

pears quite clean and well defined. However, Fig. |l^ shows 
that NGC7790 has a U-B:B-V diagram where the data 
poorly fits the ZAMS line. We show one ZAMS shifted to 
fit the OB stars which implies E(_B-V)= 0.59±0.04 and an- 
other shifted to fit the F stars which would imply E(_B- 
l/)=0.43±0.04. Most previous estimates are closer to that 
for the OB stars. This poor fit of the ZAMS to the U-B:B- 
V data in the case of this cluster is particularly significant 
since it contains three Cepheids (see Table 1). 

The U-B:B- V diagram for NGC7790 has a history of 
controversy. The original C/BF photoelectric photometry of 
Sandage (1958) of 33 ll<y<15 stars was criticised by Pe- 
dreros et al (1984). A check of 16 stars with the KPNO CCD 
seemed to confirm that Sandage's U-B and _B- F colours were 
too blue by -1-0.075 and -1-0.025 mag respectively, although 
few details were given of errors etc. However we find excel- 
lent agreement between the work here and the photometry 
of Sandage, see Table I Fig. and Fig 13. A direct com- 
parison of the photoelectric observations of Sandage to the 
photographic observations in Pedreros (see Table 1 in Pe- 
dreros et al. 1984) implies that the differences in the colours 
of Sandage are too blue by ^ 0.04 mags for both the U-B 
and B- V colours and when we compare our CCD photom- 
etry to 22 of the brightest photographically observed stars 
from Pedreros, we find similar colour differences of around 
0.04 mags in both the U-B and B- V colours. 



The U-B:B- V diagram in Pedreros et al ( 1984 ) seemed 



to the Northerly latitude of the cluster NGC129. The values 



to give the same sort of ill fitting ZAMS, throughout the 
range 0.3<-B- y<1.2, as found in this work. The suggestion 
was that the problem might lie in Sandage's photometry 
which Pedreros et al had used for calibration. However, we 
have tested various zero points for the U-B and _B- Fcolours. 
With our own zero points for the colours we find the ill fitting 
ZAMS and even if we apply the corrections suggested by 
Pedreros to the offsets found between our photometry and 
the photographic data of Pedreros we still find an ill fitting 
ZAMS. 

Romeo et al (1984) used CCD data to obtain BVRI 
photometry for this cluster to V—20. In the absence of U, 
their only route to E(_B- V) was via fitting the shape of the B- 
V:Vand V-LVCMD and they obtained E(B- VO=0.54±0.04. 
They used the Sandage (1958) photometry for calibration in 
B and V, subject to the small re calibration in these bands 
by Pedreros et al. and checked against the B, V photoelec- 
tric photometry of 10 stars by Christian et al (1985). They 
tested the faint photographic photometry of Pedreros et al 
and found scale errors at B>17 and y>15 in the sense that 
Pedreros et al were too bright. In B- V, however they claimed 
better agreement with B-Vpedreros being ~0.1mag too red. 
A comparison of our photometry and that of Romeo sug- 
gests that there is a scale error for B>17 and V>16 but 
the extent of this is less than in the comparison of the Pe- 
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Figure 6. The top panel shows the £?- y colour-magnitude dia- 
gram for the cluster NGC7790 but with only the stars that fit the 
zero age main sequence shown. The lower panel shows the same 
stars as they appear on the f7-B.B- y diagram. The U-V deficit is 
still clearly apparent among the stars that are most likely to be 
members of the cluster. 



dreros et al. data with the Romeo data, at V—17 our data 
is brighter than Romeo's by 0.1 mag whereas the Pedreros 
et al. data is brighter by 0.2 mags and similar differences 
are found in the B data. We conclude that the photome- 
try in th i s stu dy agrees very well with the photometry of 
Sandagj ( 1958) and is better agreement with the CCD data 
of Romeo et al (1989) than that of Pedreros, (see Fig. ^ for 



more details). 

We have also attempted to check to see if contamination 
IS causmg the ill fitting ZAMS. Figure | shows the B-V:V 
and the [/-5;B-1/ diagrams for the cluster NGC7790. The B- 
y.'F diagram has been trimmed so that only the stars that 
lie very close to the ZAMS remain. As these stars have the 
correct combination of distance and reddening to lie almost 
on the main sequence then it is likely that they are main 
sequence stars. The same stars are then used to produce 
the U-B:B-V colour-colour diagram. The same UV deficit 
around the F-type stars that appears in Figure when 
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Figure 7. The open circles show the photometry of this work. 
The solid triangles show the photometry of Sandage which does 
not go deep e n ough to test the ill fitting ZAMS. The results of 



Pedreros et al (1984) (filled squares), though not as deep as those 
presented here, show the same effect as the photometry in this 
work and confirms the poor fit of the ZAMS U-B:B-V re\a,tioa in 
NGC7790 



all the stars in the field are used, appears in Figure Fig- 
ure ^ shows the colour-colour diagram found from this work 
with the points from Sandage and Pedreros included. The 
Sandage points do not go deep enough to test the shape 
of the data but the Pedreros et al. points do and the poor 
match of the data to the ZAMS is seen. Therefore, we believe 
the ill fitting ZAMS to the (7-i?;B- F diagram is not caused 
by contamination from foreground or background stars or by 
errors in the photometry but is a real feature in the data. 

Thus our estimate of the NGC7790 reddening based 
on the I7-5;B- F colours of OB stars, E(B- V)=0.59±0.05, is 
between the E(B- 1/1=0.52 ±0.04 of Sandage (1958) and the 
E(B-y)=0.63±0.05 of Pedreros et al(1984) who used similar 
techniques. 

Assuming the estimate of E(B- 'l/)=0.59±0.05 from the 
OB stars, the distance modulus which fits V:B- V is then 
^o = 12.72±0.11 which is close to the original value of 
12.8±0.15 found by Sandage(1958), although this is in the 
opposite direction that would be expected from the differ- 
enc e in reddening and slightly more than 12.65 obtained 
by Romeo et al (1989) which is roughly in line with their 



obtaining E(5-l/)=0.54 for the reddening. Pedreros et al 
(1984) obtained only 12.3 for the distance modulus, no fits 
are shown in the paper and the reaso n they obtain this low 
value is unclear. The v alues quoted by Laney fc Stobie ( 1994 ) 
are similar to those of Pedreros et aj ( |1984| ). 



Sandage ( 1958| ) states that the membership of CF Cas, 
CEa Cas and CEb Cas to NGC7790 is almost certain due 
to the position of the Cepheids on the CMD. 
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Cepheid Log(P) E(B-V)ccph <V> <K> 



Mr. 



A(Mv) A{Mk) 



V367 


0.799 


1.27 


11.604 


6.662 


-3.789 


-4.937 


-3.755 


-4.992 


-0.034 


0.055 


U Sgr 


0.829 


0.45 


6.685 


3.952 


-3.841 


-5.232 


-3.781 


-5.214 


-0.060 


-0.018 


EV Set 


0.490 


0.62 


10.131 


7.028 


-2.877 


-4.164 


-2.191 


-3.553 


-0.686 


-0.611 


WZ Sgr 


1.339 


0.50 


8.023 


4.565 


-4.810 


-6.738 


-4.657 


-6.789 


-0.153 


0.051 


TW Nor 


1.033 


1.24 


11.670 


6.319 


-3.560 


-5.156 


-3.797 


-5.411 


0.237 


0.255 


QZ Nor 


0.730 


0.27tt 


8.866 


6.662 


-3.164 


-4.634 


-3.128 


-0.036 


-0.046 


-0.035 


V340 Nor 


1.053 


0.38 


8.375 


5.586 


-4.056 


-5.699 


-3.797 


-5.640 


-0.259 


-0.059 


CV Mon 


0.731 


0.84 


10.306 


6.576 


-3.821 


-5.072 


-3.305 


-4.896 


-0.516 


-0.176 


RU Set 


1.294 


0.93(1 


9.465 


5.071 


-4.901 


-6.508 


-5.186 


-6.775 


0.284 


0.267 


SV Vul 


1.654 


0.44(1 


7.243 


3.920 


-5.440 


-7.415 


-6.028 


-8.004 


0.588 


0.589 


DL Cas 


0.903* 


0.60 


8.97* 


5.93t 


-3.648 


-5.136 


-3.860* 


-5.36t 


-0.171 


0.181 


CF Cas 


0.688* 


0.55 


11.14* 


8.01t 


-3.341 


-4.872 


-3.170* 


-4.85t 


-0.171 


-0.022 


CEa Cas 


0.711* 


0.55 


10.92* 


NA 


-3.561 


NA 


-3.390* 


NA 


-0.171 


NA 


CEb Cas 


0.651* 


0.55 


10.99* 


NA 


-3.493 


NA 


-3.330* 


NA 


-0.163 


NA 



5. Comparison between the work here anc 



IhaLol 



.a.ncy gtobic ( 1994 ) . * indieates where values are taken from 



Feast & Walker 



( |1987| ) and f where the values are inferred from Welch et al (1985). H indieates where the Cepheid reddenings have been corrected (see 
text) 



4 P-L RELATION 



Using the values for the distance modulus and the redden- 
ing towards the cluster, we proceed to determine the P-L 
relation. As well as the reddening and the distance modu- 
lus of the cluster, the apparent magnitude of each of the 
Cepheids is required. Where possible these come from LS 
(1993, 1994) who have high quality V and if-band measure- 
ments for most of the Cepheids in this study. The clusters 
NGC7790 and NGC129 lie at northerly latitudes so are un- 
observable from SAAO and so there are no magnitudes from 
Laney and Stobie for these Cepheids. Th e j<"-band data fo r 
DL Cas and CF Cas comes therefore from Welch et al (11985 



and the F-band data and the periods are taken from | 
& WallQ (|1987|). 



The reddening obtai ned from the ZAMS fi tting is that 
of the cluster OB stars. Schmidt-Kaler (1982) found that 
when the effect of the colour difference between the OB stars 
and the Cepheid is taken into account, 



E(B-F)ccph = E(B-l/)ci,e[0.98- 

0.09(< Bo> - <Vo>) ccphj 



(2) 



gives a good approximation to the reddening of the Cepheid. 
The val ues fo r < Bo > — < Vo >)coph come from Feast & 
Walkei[|( |l98l ). 

Fi gure |8| shows the Mi/-Mr- - Log(P) relation with the 



reddeni pgs found for the Galactic Cepheids (filled circles) 



found from the cluster reddenings obtained in this work via 
equation ^ The triangles show the same relation for the 
LMC Cepheids and the squares are for the SMC Cepheids 



These are taken from Tables 2 and 3 in Laney & Stobie 
(1994). Three of the Cepheids, (from left to right in Fig- 
ure ^) QZ Nor, RU Set and SV Vul seem to have the wrong 
Mv-Mk colours for their periods. The clusters containing RU 
Set and SV Vul suffer both from the presence of differential 
reddening and from the fact that the Cepheids lie at some 
distance away from the cluster. Figure ^ indicates that the 
reddeni ng local to RU Set and SV Vul may be somewhat dif- 



ferent f 'om the average value of the cluster reddening. We 



correct for th is using the space reddenings given in Laney & 
StobieQ|l99|) which are more local to the Cepheids (see for 
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Figure 8. The Galactic Cepheids are shown by the solid circles, 
the LMC Cepheids by t he triangles and the S MC Cepheids by 
the squares, taken from [Laney &: Stobiej ( |l994| ). The three solid 
circles in a box are QZ Nor, RU Set and SV Vul. The re ddenings 



of these Cep heids are corrected to the values given in Laney & 
Stobie ( |l993| ), indicated here by the boxed open circles. 



example Turner 1980). Note that because there is difficulty 
obtaining the Cepheids true reddening from the cluster red- 
dening for these two Cepheids we do not include th em in 
the best sample (later in this section). Walker (1985b) notes 



that the Cepheid QZ Nor also lies away from the centre of 
the cluster NGC6067, at a distance of two cluster radii so 
again the reddening of the cluster may not be appropriate 



for the reddening of the Cepheid. Laney & Stobie ( 



199^) 



take the v alue for the Cepheid reddening from Coulson and 
Caldwel] ( |l984| ) of E(B- V) = 0.265, derived from BVIc red- 
denings. This is the reddenirig used to calculate the position 
of the open circle in Figure H and which we assume for QZ 
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Nor henceforth. However, the effect of changing the redden- 
ing of the cluster to E(5-y) = 0.265 changes the distance 
modulus to the LMC by less than 0.02, less than the error 
quoted here. 

To determine the absolute magnitude of the Cepheid, 
the apparent magnitude has to be corrected for reddening 
and distance. First of all the extinction coefflcient is re- 



quired. We follow Laney & Stobie (1994) and use 
K(ceph) = 3.07 + 0.28(B - V)a + 0.04E(B - V)coph 



(3) 



to take into account the effect of Cepheid colour on the 
ratio of total to selective extinction. The Cepheid reddening 
comes from Eq. |^ except for the three corrected values. Then 
the reddening free magnitudes of the Cepheids are 



Vo = V ■ 



Ko 



- »(ceph)E(B - y)eeph 
K(ceph)E(B 



Vo-V + K + ^ 



^)ccph 



1.1 



(4) 



The expression for Ko has the form given above as the ex- 
tinction coefficient in the Tf-band is one tenth of that in the 
V-band. To obtain finally the absolute magnitude, the dis- 
tance modulus, /io, given in Table ^ has to be subtracted 
off. 

The P-L relation can now be determined. We consider 
two samples, one where we consider all the Cepheids avail- 
able to us and another where the Cepheids RU Set and SV 
Vul are removed due to the problem of differential reddening 
and the question of cluster membership. The zero points are 
obtained by fixing the slope and obtaining the least squares 
solution using the Galactic Cepheids in this study. These are 
summarised in Table ^ The slop es that are considered are 
the slopes from Laney fc Stobiq (1994) which are the best 



fitting slopes to all the Cepheid data (Galactic open clus- 
ter Cepheids, LMC and SMC Cepheids) in their study. The 
slopes are -2.874 in the F-band and -3.443 in the K-h&nd. 
Also considered is -2.81 in the F-band as this is the slope of 



Feast & Catchpole 



t he LM C Cepheids and the slope used by 
( [1997| ). 

Once the slope and zero point of the PL relation i s fixed, 
the distance modulus to the LMC can be calculated. Laney 



Stob[e| ( |1994| ) give the period and the dereddened V and 
7f-band magnitude, Vo and Ko of 45 LMC Cepheids. The 
distance to the LMC is then given by 



< Ho >- 



IN 



N 



■ S * log(P)) 



(5) 



where ruo represents the dereddened apparent magnitude 
in each waveband and 5 and p are the values for the slope 
and zero point as given in Table ^. Our PL( and PL(it) 
relations are shown in Figure ^ alo ng with the equ ivalent 
relation for the LMC Cepheids from Laney & Stobie (1994) 
using the distance moduli in Table ^ to determine the abso- 
lute magnitude. 

Taking the value for the PL( V) zeropoint for the best 
sample with the -2.874 slope used by Laney & Stobie (1994) 
gives p=-1.279±0.33 which is in good agreement with the 
value of p=-1.197±0.09 found by these authors and implies 
a distance modulus of 18.57 for the LMC as compared to 
18.50 found by Laney & Stobie. Figs. 8(a,b) show that the 
best fitting zeropoint gives a somewhat poor fit to the ma- 
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Figure 9. The Cepheid P-L relation, (a) shows the F-band P- 
L relation with a slope of -2.874 and the (b) shows the F-band 
P-L relation with a slope of -2.810. (c) shows the K-band P-L 
relation with a slope of -3.443. The solid symbols are the Galactic 
Cepheids, marked specifically are U Sgr (diamond), EV Set (filled 
triangle), CF Cas, CEa Cas and CEb Cas (filled squares). There 
is no iiT-band data for the Cepheids CEa Cas and CEb Cas. The 
zeropoints and slopes are shown in each panel. The open triangles 
show the P-L relation for the LMC Cepheids. The periods and 
magn itudes are taken from Table 2 and 3 in Laney & Stobie 



(1994) and the distance moduli to the LMC are the 'best' values 
from Table M in this work. 
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Band 


Sample 


Slope(<5) 


Zeropoint(p) 


Ato (LMC) 


y 


All(14) 


-2.874 


_i 22q-l-n 34 




V 


All(14) 


-2.810 


-1.289±0.33 


18.50±0.036 


K 


All(12) 


-3.443 


-2.148±0.30 


18.42±0.019 


V 


Best(12) 


-2.874 


-1.279±0.33 


18.58±0.036 


V 


Best(12) 


-2.810 


-1.332±0.32 


18.55±0.036 


K 


Best(lO) 


-3.443 


-2.200±0.29 


18.47±0.019 


V 


logP<1.0(9) 


-2.874 


-1.418±0.19 


18.71±0.036 


V 


logP<1.0(9) 


-2.810 


-1.465±0.19 


18.68±0.036 


K 


logP<1.0(7) 


-3.443 


-2.314±0.21 


18.58±0.019 



Table 6. The zeropoints for the P-L relation and distance mod- 
ulus to the LMC 



jority of the Galactic Cepheids which lie at log(P)<1.0. This 
is because the 12 Galactic Cepheids in the best sample give 
a slope of 5=-1.85±0.33 which is much flatter than the LMC 
data which gives (5=-2.79±0.1, close to -2.81. Indeed, if only 
the 9 Galactic Cepheids with log(P)<1.0 are used, the zero- 
point rises to p=-1.418±0.19 and the LMC distance modulus 
would rise to 18.70, indicating why the errors on the PL( V) 
zeropoint are as large as they appear in Table 5. 

Our PL( V) zeropoint is also consistent with the ze- 
ropoint and LMC distance obtained from an analysis of 
Hipparcos tri gonometrical paralla x es of nearby Galactic 
Cepheids by Feast & Catchpole (L997). They obtained 
p=-1.43±0.1 for the Galactic PL( 1/) zeropoint for an as- 



sumed slope of (S=-2.81. This can be compared to the p—- 
1.332±0.32 obtained for our best sample with the same 
slope. They used the same 45 Laney & Stobie (1994) 
Cepheids as used here to obtain a metallicity corrected LMC 
distance modulus /io=18.70±0.10. We note that their semi- 
theoretical metallicity correction to the LMC Cepheid V 
magnitudes increases the distance to the LMC, which is in 
the opposite sense to most empirically determined estimates 
of the efTects of metallicity on Cepheids (eg Kennicutt et 
al 1998). Subtracting their metallicity correction leads to 
an LMC distance modulus /^o = 18.66±0.10 which can be 
directly comnared with our best value of /io=18.55±0.036 
from Table m. We conclude that our PL( V) estimates the 
LMC distance modulus are between those of Laney & Sto- 
bie and Feast & Catchpole but have too little statistical 
power to discriminate between these previous estimates. 

The K-haxiA P-L relation is tighter for the LMC 
Cepheids and for the Galactic Cepheids, see Figure ^(c), 
and the slopes are closer with the LMC Cepheids giving 5=- 
3.27±0.04 and the best sample of Galactic Cepheids giv- 
ing ^=-2.81±0.21. We note in passing that the slope of the 
Galactic Cepheid PL(7t) relation is now much flatter than 
the 5=-3.79±0.1 slope found in the Galactic Cepheid sam- 
ple of Laney & Stobie. Assuming the -3.443 slope used by 
Laney & Stobie our best sample in Table 5 gives a PL(it) 
zeropoint of p=-2.200±0.29 which implies an LMC distance 
of /io=18.47±0.29 which remains in good agreement with 
the value Mo=18.56±0.07 found by Laney & Stobie. (1994). 
The smaller error of Laney & Stobie is due to their larger 
numbers of calibrators although it must be said that many 
of their extra calibrators (8/12) are in associations rather 
than clusters and frequently given half-weight in P-L fits. 
Indeed, Hipparcos proper motion data has shown that one 
of their further cluster Cepheids, S Nor, is also unlikely to 
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Figure 10. The approximate positions of all the clusters con- 
sidered here. The clusters NGC6664 (triangle), M25 (diamond), 
NGC129 and NGC7790 (square) are highlighted. These symbols 
are the same as in Figure tJ 



be a member of its cluster (Haguenau conference, Septem- 
ber 1998). Moreover, they have not included the 4 cluster 
Cepheids in NGC129 and NGC7790. Therefore we believe 
that our result supercedes the Laney & Stobie result with 
our bigger error estimate perhaps being a more realistic in- 
dication of the actual errors. Certainly in our best sample 
the biggest changes in Mx between Laney & Stobie and our- 
selves, which contribute most to our 50% increased scatter, 
are for EV Set and TW Nor (see Table 4) where the reasons 
for the distances used by Laney & Stobie are unclear. 

Finally, as our overall estimate of the LMC distance, 
we take the average of the PL( V) and PL(i(^ estimates in 
the best sample of Table 5 which gives /io = 18.51±0.3. We 
conclude that although in the case of individual clusters we 
have markedly improved the distance and reddening esti- 
mates, our new estimates of the zeropoint of the PL relation 
and thus the distance to the LMC are close to previous val- 



5 DISCUSSION 

We now discuss the most intriguing new result in this study, 
which is that the Solar metallicity ZAMS may not always 
fit the (7-5:5- F data in individual clusters. This is not the 



first time an effect like this has been seen. Turnei (1993) saw 
poorly fitting U-B:B-V TjAMS for the open cluster Roslund 
3. Turner interpreted this as evidence for the young, B-type 
stars having a cocoon of circumstellar dust around them. 
This cocoon of dust then increases the reddening of the B- 
type stars as compared to the F and G type stars, causing 
the ill-fitting [/-5;5- VZAMS. The O and B type stars in the 
i7-5.'5-F diagram of Roslund 3 show a large spread around 
the ZAMS as the amount of excess dust would probably 
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vary from star to star. Excess reddening could perhaps also 
be so strong that it was causing some O and B stars to be so 
reddened that they appeared as F type stars. However, the O 
and B type stars in the [/--B/B- V diagram for NGC7790 are 
very tight so the shape of the U-B:B-V diagram is unlikely 
to be caused by excess dust. 

The next possibility we consider is that the effect might 
be due to stellar evolution, However, the CMD for the clus- 
ters look unevolved even at AOV as might be expected for 
clusters which have Cepheid variables and are expected to 
be less than 10* years old. 

We also considered whether the discrepancy between 
the main sequence fitted distance and the Hipparcos parallax 
to the Pleiades could explain our result. If it were assumed 
that all open clusters had roughly the same composition then 
the different forms for U-B:B-V that we find might be taken 
as evidence that the colours of main sequence stars may not 
be unique. This possibility has also been discussed as an ex- 
planation of_th£_groblem_with the MS fitted distance to the 
Pleiade s (|van Leeuwen 1999|) and if it proves relevant in that 



case it vill certainly also be worthy of further consideration 
here. 

If the reddening vector in U-B:B-V varied as a function 
of Galactic position then this would also affect our results. 
However, at least in the case of NGC7790 it seems that for 
whatever relative shift in U-B and B-V, the ZAMS still has 
the wrong shape to fit the observed colour-colour relation. 

The final possibility is that metallicity is affecting the 
F stars' U-B colours in some of these clusters. Qualitatively 
there is some evidence supporting this suggestion. First, 'line 
blanketing' is well known to redden the U-B colours of metal 
rich stars at F and G and low metallicity sub-dwarfs are 
known to show UV- excess as the reverse of this case (e.g. 
Cameron (1984)). Second, there is some suggestion that the 
cluster, NGC7790, that shows a UV excess lies outside the 
solar radius while NGC6664 which is redder in U-B at F 
tends to lie inside (see Figure |l^) . M25 which lies closest to 
the Sun also fits the solar metallicity f/-_B:_B- F diagram as 
well as any of the clusters. Given that metallicity in the 
Galaxy is known to decrease with Galactocentric radius, 
this is suggestive of a metallicity explanation. Many of the 
other clusters' C/-_B.-_B-y diagrams are either too noisy due to 
differential reddening (Tr35, NGC6823) or too obscured to 
reach the F stars (NGC6649, Lynga 6, vdBerghl) to further 
test this hypothesis. However, NGC6067 forms a counter- 
example to any simple gradient explanation, since it seems 
to have a normal UBV plot and lies inside the solar position 
This would have to be accommodated by allowing a sub- 
stantial variation on top of any average metallicity gradient. 

However, quantitatively the case for metallicity is less 
clear. The size of the UV excess seen is much larger in the 
case of NGC7790 than expected on the basis of previous 
metallicity estimates of these clusters, or of any measure- 
ment of the amplitude of the Galactic me tallicity gradi - 
en t. Using the Fe/H vs A U-B relations of Carnev J 197£ I 



or |Caneron| (|l985|) it would be concluded that NGC7790 
showed A U-B ~ 0.2mag which corresponds to Fe/H~-1.5. 
Thus clusters which on the basis of their Main Sequences 
and the presence of Cepheids, must be less than 10* yr 
old, wo uld be implied to h ave n ear halo metallicity. Pre- 
viously, Panagia and Tossi (1991 ) fin d Fe/H~-0.3 for these 
2 clusters. Also |Fry fc Carney| ( [19971) find Fe/H=-0.2±0.02 



for NGC7790 while finding Fe/H=-0.37±0.03 for NGC6664 
based on spectroscopy of the Cepheids in these clusters 
themselves. Also according to the Galactocentric metallicity 
gradient whic h is usually taken to lie in the range -0.02- 
O.ldex kpc~^ (Rana 1991), there should only be on the av- 
erage A Fe/H ~ 0.3 in the range of metallicity covering these 
clusters. 

On the other hand, it should be noted that Panagia & 
Tosi's Fe/H estimates are based on more poorly measured 
estimates of UV excess than those presented here and also 
that there is little agreement between the metallicity es- 
timates of Fry and Carney and those of Panagia & Tosi. 
Measuring the metallicity of the Cepheids themselves as at- 
tempted by Fry and Carney is difficult since the effective 
temperature is a function of the light curve phase and a 
small difference in estimated temperature can make a large 
difference in metallicity. Also the Galactocentric metallic- 
ity gradient at least as measured for open clusters depends 
on relative ly poo r U-B p hotometry at th e limi t of previous 



data from Janes (197S) 



Cameron 



(1985) and Panagia and 



Toss: (1991). In any case, it is well accepted that the dis- 



persion in metallicity around the mean gradient is indeed 
high, wi th the r a nge - 0.6<Fe/H<-|-0.3 at the Solar position. 
Further [Geislei] ( |l987| ) using Washington photometry to es- 
timate metallicity also found an example of a cluster, NGC 
2112, only ~ O.Skpc outside the solar radius with Fe/H=-1.2, 
although this cluster is older than those discussed here. 

However, it would also seem that the tightness of the 
P-L relations in Fig. ^ could form a final argument against 
the idea that NGC7790 has Fe/H~-1.5. If the metallicity of 
the cluster NGC7790 was really Fe/H~-1.5 then at given 
B-V, MS stars would be sub-dwarfs with ~lmag fainter ab- 
solute V magnitudes than normal solar metallicity main se- 
quence stars (see Cameron, 1984, Figs. 5,6). Thus since we 
have used a normal Main Sequence to derive the distance 
to NGC7790, is it not surprising that the Cepheids in these 
clusters lie so tight on the P-L relation when they should 
be a magnitude too bright if the low metallicity hypothesis 
is correct. The only way that the low metallicity hypothesis 
for NGC7790 could survive this argument is if it were pos- 
tulated that the effect of metallicity on Main sequence star 
magnitude and Cepheid magnitude were the same - then the 
effect of our derived distance modulus being ^1 magnitude 
too high would be cancelled out by the fact that the Cepheid is 
actually sub-luminous by 1 magnitude because of metallicity 
which would leave the Cepheid tight on the P-L relation as 
observed. This might not be too contrived if a low metallicity 
Cepheid prefers to oscillate about its subdwarf, rather than 
solar metallicity, zero-age luminosity (at fixed effective tem- 
perature) on the Main Sequence. This would lead to a strong 
implied metallicity effect on the Cepheid PL( V) and PL (if) 
zeropoints; the implication would be that 



Jf^ -0.66 in 
the sense that lower metallicity Cepheids are fainter. This 
coefficient is within the range that has been dis cussed for the 
em pirica l effec ts of metallic ity on Cepheids by Kennicutt et 
a] (|l99j) and [Gould (|l994| ) although the most recent work 
by [Kennicutt et al[ ( 199^ appears to give a lower value of 
sfT/h — 0-24zfc0.16 again in the same sense. 

The immediate effect on the distance to the LMC with 
Fe/H=-0.3 is that our estimate of its distance modulus 
would decrease from 18.5 to 18.3. However, since all that 
is determined at the LMC is the slope of the P-L relation. 
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then the zeropoints we have derived in Table 5 from the 
Galactic Cepheids would still refer the P-L relation to the 
Galactic zeropoint. The ultimate effect on Ho would then 
be decided by the metallicity of the Cepheids, for example, 
in the galaxi es observed by the HS T for the Distanc e Scale 
Key project Ferrarese et al. (199£), by Fanvi r et a l (1991:) 

Saha| ( |199q ) in the 



in the case of the Leo I Group and by 



case of SNIa. Zaritsky (1994) have measured metallicities 



in these galaxies already but if the dispersion in Cepheid 
metallicity is as large as it is implied to be in the Galaxy 
then there may be some signature in a wider dispersion in 
the Cepheid P-L relations in at least the high metallicity 
cases. The possibility of detecting this signature is currently 
being investigated (Shanks et al 2000 in prep.) 

Obviously the most direct route to checking the metal- 
licity explanation for the anomalous behaviour seen in the U- 
B:B- V diagrams is to obtain medium-high dispersion spec- 
troscopy for a sample of F stars in NGC7790 and NGC6664 
to determine the metallicity directly for these main sequence 
stars. Currently proposals are in to use WHT ISIS spectro- 
graph for this purpose. 



cluster to cluster may be affecting the U-B colours of F- 
and G-type stars. The problem is that the metallicity vari- 
ations this would require are much larger than expected for 
young, open clusters. More work is therefore required to de- 
termine the metallicity of the individual main sequence stars 
in each of the clusters NGC7790 and NGC6664. If metalhc- 
ity is proven to be the cause of the anomalous U-B:B-V 
relations, then it would imply that the Cepheid P-L relation 
in both the visible and the near-infrared is strongly affected 
by metallicity. 
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6 CONCLUSIONS 

We have presented colour-colour diagrams and colour- 
magnitude diagrams of a sample of galactic clusters which 
contain or are associated with Cepheids. All the clusters 
have been observed using similar methods and the data re- 
duction and extraction has also been done with similar tech- 
niques. The use of the improved fZ-band data has allowed 
powerful new checks of previous E(_B-V) estimates over a 
wide range of magnitudes. In order to estimate the redden- 
ings and distance moduli, we have fitted all the clusters in 
the same manner and have not attempted to correct for dif- 
ferential reddening but instead taken a simpler approach 
and fitted the average reddening value of the cluster. In 
most cases the differences that we have found between val- 
ues for the reddening and distance modulus are small and 
where there are significant differences these can mostly be 
explained by comparing whether the ZAMS fit was made to 
the centre or to the edge of the colour-colour and colour- 
magnitude diagrams. 

The Cepheid P-L relations found from fitting the best 
sample are Mv=-2.81xlog(P)-1.332 and Mx=-3.44xlog(P)- 
2.20 and a distance modulus to the LMC of 18.54±0.32 in 
the F-band and 18.46±0.29 in the i("-band giving an overall 
distance modulus to the LMC of 18.50±0.3, ignoring any 
possible effect of metallicity. These results for both the PL 
relations and the LMC distance are consistent with the pre- 
vious results of Laney & Stobie (1994) although the im- 
proved distances and reddenings have increased the errors 
over what was previously claimed. These increased errors 
mean that our result for the PL( V) relation are also consis- 
tent with the result of Feast & Catchpole (1997) from Hip- 
parcos measurements of Cepheid parallaxes, although this 
gives rise to an LMC distance modulus of fj,o= 18.66±0.1 as 
opposed to our /io=18.51±0.3. 

With the improved [/-band data, we find that for at 
least two of the clusters, the data in the [/-B.-B- F two-colour 
diagram is not well fitted by the solar metallicity ZAMS. 
One possibility is that significant metallicity variations from 
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Figure 11. The !7-S;B- F diagrams for tlie clusters used in the study. The arrow indicates the direction of the reddening vector. 
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Figure 12. The ^--B; V diagrams for the clusters used in the study 
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Figure 13. The F-Jf; F diagrams for the clusters used in the study 
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Figure 14. A comparison between tlie zero point compared in this work with the zero point obtained in previous photoelectric studies. 
The full reference for the comparison is given in Table M. Only the clusters where an independent zero points was obtained are shown. 
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Figure 15. A check for any colour dependent relationship between the work in this study and previous photoelectric work. The full 
reference for the comparison is given in Table w. Only the clusters where an independent zero points was obtained are shown. 
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Figure 16. The zero points and check against colour for the clusters where previous work had to be relied upon for calibration purposes. 
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Fig ure 17 . Comparison of all the different sources of photometry for NGC7790. 
et al | (Hs^ ) 



Note there is no U-hand data in the study by Romeo 
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